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The thermopower of superconductors measured via the magnetic flux in bimetallic loop is evaluated. 
It is shown that by a standard matching of the electrostatic potential, known as the Bernoulli 
potential, one explains the experimentally observed amplitude and the divergence in the vicinity of 
the critical temperature. 



The thermopower is a widely used tool to study elec- 
tronic properties of conductive materials. An exception 
are superconductors, where the supercurrent cancels any 
diffusive current so that the zero net current or voltage 
are observed. This feature is known from 1935 [1] and 
since then there were a number of attempts to access 
diffusive currents in an indirect way. 

Already in 1944 Ginzburg [2] noticed that in inho- 
mogeneous systems like bimetallic loops, the counter- 
flowing supercurrent creates a magnetic flux. The boom 
in this field came thirty years later. In 1974 Garland and 
Van Harlingen proposed a simple phenomenological the- 
ory [3] and Gal'perin, Gurevich and Kozub published a 
microscopic treatment [4] based on the Boltzmann-type 
approach. These theories predicted fluxes of similar am- 
plitudes and temperature dependences. 

In the same year Zavaritskii presented first experimen- 
tal data [5] and he was soon followed by others [6-8]. 
Experimental results were a surprise. Zavaritskii [5] and 
Falco [7] observed the expected temperature dependence, 
but Pegrum, Guenault and Pickett [6] and Van Har- 
lingen and Garland [8] monitored a thermally induced 
magnetic flux by five orders of magnitude larger. More- 
over, the theory predicts that close to T c the flux $ di- 
verges as oc (T G — T) _1 , while a steeper divergence 
oc (T c — T)~ 3 / 2 was observed [6,8]. The experimental 
situation in the late 1970 is reviewed in Ref. [9]. 

The giant flux stimulated a number of theoretical stud- 
ies [10-15] that explored various additional components 
ranging from a trapped flux, over impurities, over in- 
terfaces, to an influence of supercurrent flow. Most of 
these ingredients bring only a minor correction to the 
original prediction. It was speculated, that the only siz- 
able contribution can come from the trapped flux, which 
increasingly leaks into the ring as the temperature ap- 
proaches its critical value. All these speculations were 
terminated by measurements of Van Harlingen, Heidel 
and Garland [16]. To avoid the penetration of the exter- 
nal magnetic field they used the toroidal geometry and 
convincingly demonstrated that the large magnetic flux 
with the ^ oc (T c — T)~ 3 / 2 divergence is a genuine effect. 
By comparing a number of samples they could conclude 
that the flux is proportional to the thermopower in the 
normal state and therefore that it is indeed caused by the 



thermal diffusion of electrons. 

The lack of at least a qualitative theory has discour- 
aged further measurements in this direction and the ther- 
mopower joined the family of puzzling transport proper- 
ties in superconductors. Alternative measurements of the 
thermopower via the superconducting fountain effect [17] 
or the charge imbalance in the conversion region [18] yield 
theoretically expected values but they have a rather low 
accuracy. In result, the values of the thermoelectric co- 
efficients of superconductors are still not available, what 
contrasts with extensive data gathered for superconduct- 
ing materials above T c . 

New theoretical interest in the giant magnetic flux 
emerged after seventeen years. Marinescu and Over- 
hauser [19] have analyzed the theory of Gal'perin, Gure- 
vich and Kozub [4] and concluded, that its failure indi- 
cates a conceptual mistake in the underlying Boltzmann 
type transport theory developed by Bardeen, Rickayzen 
and Tewordt [20]. They made an ad hoc modification 
of the transport theory by including the momentum ex- 
change between the condensate and quasiparticles. With 
this modification, a good agreement between theory and 
experimental data was reached. 

The modified transport theory, however, is in conflict 
with other properties of superconductors as discussed re- 
cently by Gal'perin, Gurevich, Kozub and Shelankov [21]. 
From the time-reversal symmetry they showed that this 
theory predicts dissipative currents also in equilibrium 
systems with inhomogeneous chemical composition, i.e., 
in any real superconductor. 

Here we show that explanation of the thermopower 
does not require any changes in the theory of transport 
in superconductors. The legitimate request of Marinescu 
and Overhauser to cover properly the balance of forces 
between the superconducting and normal electrons is nat- 
urally covered by the theory of electrostatic potential in 
superconductors known as the Bernoulli potential. Un- 
like forces ad hoc added to the kinetic equation, forces 
derived from a scalar potential cannot result in an artifi- 
cial dissipation. 

Below we demonstrate that the giant magnetic flux 
can be described in a simple manner with the help of the 
Bernoulli potential. Our approach parallels the textbook 
theory of thermopower in normal metals in that we eval- 
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uate the net current in the sample from the requirement 
of the electrostatic potential matching. 

We will use notation related to the experimental setup 
of Van Harlingen, Heidel and Garland [16] shown in 
Fig. 1. The sample is a toroid with the internal cylin- 
der from the Lead and the external one from the Indium. 
The magnetic flux $ in question is restricted to the vol- 
ume between cylinders, i.e., it is inside the toroidal cavity. 
The diamagnctic current J corresponding to the flux <f> 
flows on the outer surface of Lead and the inner surface 
of Indium. The flux is linear in the current, 



$ = J —fi Dln , 

2ir r Pb 



(1) 



where D is the length of the sample, and n n ,pb are the 
radii which enclose the flux. 
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FIG. 1. The toroidal sample of Van Harlingen at al [16]. 
The measured flux <E> is in the toroidal cavity. 

The experiment is aimed to measure the transport co- 
efficient Lt which determines the diffusive electric cur- 
rent jdif = — LtVT caused by the temperature gradient 
VT. As already mentioned, jdif is not observable since it 
is cancelled by the supercurrent j = —jdif- The theory of 
Garland and Van Harlingen [3] uses the London gauge to 
find the vector potential A = — ^ A 2 j, where A is the Lon- 
don penetration depth. The flux is then an integral along 
the bimetallic loop, $g = § drA = —fj, § dv\ 2 L T VT. 

The London gauge is justified only for small fluxes. 
The data from [16] show, however, that the flux is large 
so that it is of form $ = + $q , where N is an integer 
quantum number of the superconducting condensate and 
<&o is the elementary flux. Estimates [16] indicate that 
$G -C $0, therefore to understand the giant flux, we have 
to find which state N is the most favorable for the system 
with the imposed temperature gradient. Measured values 
of N in [16] range up to 250. These values are sufficiently 
large for the classical approximation, where N is treated 
as a continuous variable. Accordingly, we will not assume 
quantum restrictions of the flux <F 

The mechanism by which the flux arises is as follows. 
The diffusive current generates magnetic field, which is 
screened by the counterfl owing supercurrent. In the sur- 



face layer of the London penetration depth thickness, the 
cancellation is not complete. Accordingly, the supercur- 
rent density is a sum of the counter-flow —jdif and a 
missing counter-flow je~ x ^ x , where x is a distance from 
the surface enclosing the cavity. We call it a diamagnetic 
current, as it screens the bulk of superconductor from the 
magnetic field, which is present in the toroidal cavity. 

Our aim is to find amplitudes j in the Indium and 
the Lead. The total current J — J cPr(j + jdif) is the 
integral of the current density across each cylinder, i.e., 
J = 27rri n Ai n ji n , and due to continuity condition also 
J = — 27rrpbApbjpb- Since A depends on the tempera- 
ture, the surface values of the diamagnetic current densi- 
ties ji n ,Pb change along the temperature gradient, while 
the product Xj stays constant. 

Now we specify the condition for the total current 
J from the requirement of the scalar potential match- 
ing. As was observed by Bok and Klein [23] and with a 
higher precision by Morris and Brown [24,25], current in 
the superconductor induces perpendicular electric field. 
It is well approximated by the electrostatic potential of 
Bernoulli type [22] 



n s 1 2 
n 2 



(2) 



where n s is the density of superconducting electrons. The 
velocity of the superconducting electrons at the surface 
is given by the current density, 



v = 



L T VT 
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where plus applies for the Indium and minus for the Lead, 
in which the total current flows in opposite direction. The 
first term is due to the compensating supercurrent — jdiff , 
the second term is caused by the diamagnetic current j. 

The electrostatic potential has to be continuous, there- 
fore the potential differences created by the temperature 
gradient in the Lead and the Indium has to be equal, 



fPb(T 2 ) - </?Pb(Ti) = pin(T 2 ) - (MTl). 



(4) 



This is the central equation in our approach. From the 
set (2-4) one can directly evaluate the current J and the 
magnetic flux (1). 

Condition (4) is a simple quadratic equation, which 
includes material parameters of both, the Lead and In- 
dium. The experiment [16] explores temperatures close 
to critical temperature of Indium T^ n = 3.4 K, therefore 
only a small fraction of electrons remain superconducting 
nl n -c n In in the Indium arm. Since the critical tempera- 
ture of Lead Tjf b = 7.19 K is considerably higher, the ma- 



jority of electrons are superconducting, n 
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and 



consequently the difference of Bernoulli potential in the 
Lead is much smaller than the potential difference in the 
Indium. Briefly, the Lead effectively short-circuits ends 



2 



of the Indium, so that (4) reduces to </?in(T 2 ) = <pi n (Ti) 
and material parameters of Lead drop out. 

The second simplification follows from the relation 
between the superconducting density and the London 
penetration depth, A 2 = m/(/ioe 2 n s ). This allows 
us to express the condition (4) on the potential as 
win(Ti)/Ai n (Ti) = -win(T , 2)/Ai n (r 2 ). The net current 
is now a solution of a linear relation and the resulting 
magnetic flux reads 



1 



H L T (T 2 - Ti) (A(Ti) + A(T 2 ))n n ln 



(5) 



In this expression, the thermoelectric coefficient Lt and 
the London penetration depth A are of the Indium. To 
make the expression compact we write DVT = T 2 — T\ . 
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FIG. 2. Thermally induced magnetic flux in a toroid. The 
dots are the experimental data of Van Harlingen, Heidel and 
Garland [16] , the line is according to formula (5) of this paper 
with Ti = T c - 7mK and T 2 = T. 

In Fig. 2 we compare experimental data of Van Har- 
lingen, Heidel and Garland [16] with formula (5). The 
diameters of toroid are n n = 3 mm and n n = 1 mm. 
Material parameters are the thermoelectric coefficient of 
the normal metal, Lt — 9.82 10 3 A/Km, and the London 
penetration depth at the zero temperature, Ao = 400 A. 
All the parameters are from Ref. [16] as values for the 



sample T-4. The only open question is the temperature 
dependence of the London penetration depth close to the 
critical temperature T c In = 3.4 K. We use the asymptotic 
BCS relation A = A / ^2 - 2T/T C . 

In the linear region, T\ — T 2 <C T c — T 2 , the theory 
agrees with the experimental data within experimental 
errors. Perhaps we should return to the original aim of 
the measurement and conclude that experimental data 
from [16] confirm that the thermoelectric coefficient Lt 
close below T c has the same value as close above. 

In the non-linear region the theory deviates from data. 
This is no surprise since the presented theory is locally 
linear in the temperature gradient. Moreover, additional 
non-linear effects are caused by the so called thermo- 
dynamic correction to the electrostatic potential [26,27]. 
We aim to discuss these corrections in a next paper. 

We should mention that formula (5) has been de- 
rived under a tacit assumption that a width wi n of the 
Indium cylinder is sufficiently larger than the London 
penetration depth A. For A(T 2 ) — > Wi n the magnetic 
flux given by formula (5) approaches $ — > 2 $n, where 
$ n = -fioLrffi — Ti)wi n ri n In The flux $ n develops 
when the Indium makes transition to the normal (non- 
superconductive) state, while the Lead remains super- 
conducting. The factor i results from the unrestricted 
integration of the diamagnetic current into the bulk of 
Indium valid only for A <C w. Assuming the upper in- 
tegration limit, J = f™ dxj, one finds that $ — * $ n . 
Unfortunately, details of the flux in the narrow vicinity 
of T c have not been measured. One can merely speculate 
that <E> n is the upper limit of the diverging flux <!>. 

Finally, we want to clarify the simple potential match- 
ing used above. First, the potential has to match across 
the whole sample while equation (4) was obtained by 
matching only at the inner surface. At the outer surface 
the Bernoulli potential (2) is zero everywhere so that 
the matching is clearly satisfied. The potential profile 
between the inner and outer surfaces is nontrivial since 
the current profile is complicated by itself. Indeed, the 
screening currents in the Indium and in the Lead spread 
over different London penetration depths, and they have 
to match across the whole interface with the current con- 
tinuity satisfied at each point. We plan to evaluate the 
current and potential profiles in future. The present the- 
ory is based on our believe that the matching at the sur- 
face points of the interface is sufficient. 

Second, we have ignored the role of the flat pieces in 
the upper and lower end of sample. In these pieces, the 
temperature gradient is absent, nevertheless, the poten- 
tial difference across each piece is nonzero since the cur- 
rent density at the matching corner to the Lead is higher 
than the current density at the Indium corner. Sending 
VT to zero in (3) and using obtained velocities in (2), we 
find that the potential differencies across the upper and 
the lower flat pieces are identical and thus cancel. 



3 



Third, the Bernoulli potential (2) is the simplest ap- 
proximation of the electrostatic potential. Why we ignore 
more sophisticated potentials that include the thermody- 
namic corrections [26] and non-local corrections due to 
the finite Ginzburg-Landau coherence length [28]? Both 
these corrections result in a surface dipolc [29] which 
makes the potential matching more complex. On the 
other hand, with the surface dipole there also appears a 
dipole at the interface of Indium and Lead. We expect 
that these dipoles tend to cancel in the final matching 
condition. 

In conclusion, we would like to encourage measure- 
ments of thermoelectric coefficients in superconductors. 
Since the detection of magnetic fluxes is extremely sensi- 
tive and fluxes 10~ 3 $o can be conveniently monitored, it 
should be possible to access Lt in a wider temperature 
region, not merely few mK's bellow T c . 

Far from T c the present theory is not valid, since the 
magnetic flux becomes small and it has to exhibit the 
quantization. Fluxes smaller than the elementary flux 
are covered by the former theory [3,4], as confirmed by 
Zavaritskii [5] and Falco [7] , who monitored fluxes of the 
order of 10 -2 <J>o and 10 _1 <J>o, respectively. 

Interesting features might appear at the intermedi- 
ate region, where fluxes are comparable to the elemen- 
tary flux. For instance, the above discussed sample has 
the classical estimate of thermally induced flux equal to 
l<I>o/mK at the temperature T = T c — 65 mK. It should 
be thus in an access of experiment to observe whether 
the flux increases in steps $o or smoothly. 

We are grateful to Van Harlingen for a kind permis- 
sion to reproduce their figures. This work was sup- 
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